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Abstract 

In this report I attempt to outline the process of developing and build- 
ing an absorption imaging system capable of imaging ultracold atoms. In 
the theory section I will discuss the elements required to estimate the atom 
number and derive the equations required to determine the temperature of 
the atom cloud. I will outline the general features of a microscopic - two 
lens imaging setup via the principles of geometric optics. A discussion of 
optical design process and its features will be carried out with an example 
on absorption imaging. The essence of basic optical design process will be 
given. 
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1 Introduction 

An optical imaging system is at the heart of almost every physics experiment. 
Efficient collection of light by an optical system and its transformation into an 
image can provide a great insight into various physical processes that take place. 

Discovery and advancement in laser cooling techniques has led to a rapid de- 
velopment of ultra - low temperature physics where new states of matter were 
observed for the first time since it's prediction by Bose and Einstein in the early 
20th century. 

One of the most striking predictions of quantum mechanics and Bose - Einstein 
statistics is the phenomenon of Bose - Einstein condensation, where a dilute and 
weakly interacting gas confined in an external potential would display quantum 
effects on a macroscopic scale. The need to observe such phenomena and test 
these theories has led to the development of various optical imaging techniques 
ranging from absorption to fluorescence and more recently light sheet detection. 
In this report however, we will only discuss absorption imaging. 

Absorption imaging is by far the most popular imaging technique used in ultra- 
cold atom experiments. It can provide information about the spatial distribution 
of atoms, the atom number and its average temperature. [1] 



2 Theory 

2.1 Absorption imaging and estimation of atom number 



In absorption imaging a cold atom cloud is illuminated by a laser beam which 
is resonant with it. The absorption of light by the atoms casts a shadow which 
is then imaged onto a CCD camera. The amount of light absorbed by the cloud 
is the optical density and the corresponding effect on the intensity of the laser 
beam is given by Beers Law: 



Where, Jo, is the intensity of the incident laser beam and, J, is the intensity 
of outgoing laser beam. If the incoming intensity is below saturation, then we 



I = I e 



-OD 



(1) 
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can rewrite equation (1) in terms of the atom column density, n(x,y), and the 
absorption cross section, a [1]: 

/ = I e- n(x ^ (2) 
Here the absorption cross section is given by: 

Where, <7 , is the on resonance cross section and is given by [2]: 

Here, A is the detuning of the laser field from the atomic resonance, T, is the 
natural decay rate of the excited state and, I sat , is the saturation intensity. 

Now we can equate equations (1) and (2): 

OD(x,y) = n(x,y)a (5) 

Experimentally, we can measure the optical density by measuring the incoming 
and the outgoing intensities of the laser beams. However, to measure the actual 
optical density, we need to take into account the effects of saturation and dark 
currents as well as stray light. Taking all this into account, the final optical 
density is calculated using [3]: 

OD actual = OD mod + (1 - e ° D ^)^ ( 6 ) 



With OD mod given by: 



1 _ e -OD sat 

0D ™ d = ln e -OD meas _ e - ODsat ( 7 ) 



And finally, the OD meas is given by: 



flight Idark 



OD meas = In ; n9ni _ p (8) 

-L atoms ^dark 

Where, Iu g ht, is the intensity of light without the atoms, Idark, is the intensity 
of the probe beam off and camera fully open and, I atoms, is the intensity of light 
with atoms. 

Now, to estimate the atom number we rearrange the equation (5) and integrate 
along the x and y directions to get the total atom number: 

/ n(x,y)dxdy = - / OD actua idxdy (9) 

-oo J~oo ® J —CO J —oo 

Here the integration of the optical density (OD) implies that it's intensity com- 
ponents, J, are integrated across the image since the intensity varies spatially. 
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2.2 Derivation of the temperature of the cloud 

We assume a non-interacting cloud that obeys Boltzmann's statistics. Such par- 
ticles thus follow the Boltzmann distribution which states that the probability of 
a particle being in the Ef h energy state at temperature, T is given by: 

1 Ei 

Pi = ^ ex P(--j^f) dE ( 10 ) 
Where, Z, is the partition function and, ks, is the Boltzmann constant. 

Now suppose we have N such particles in the cloud, then, the mean number 
of particles in the E\ h energy state at temperature, T is given by: 

N E 

n t = —exp{--^-)dE (11) 
Z kb± 

We assume that these particles have motions that can fit the classical picture. 
Thus, we can express the kinetic energy of such particle as: 

E k = l -mv 2 (12) 

We would now like to express our number distribution in terms of velocity differ- 
ential instead of energy, this can be achieved by the following trick: 

1 dE 

Ek = -mv 2 =>- — — = mv =>■ dE = mvdv (13) 

2 dv 

Thus, by replacing E and dE in terms of velocity components and differentials 
we obtain: 

NmV mv 2 

rti = — - — exp{ — —)dv (14) 

Z FK 2k B T J v ; 

In the time of flight experiments, (TOF), we cannot directly measure the velocities 
of each individuals atoms in the cloud. However, we can measure the overall 
expansion of the cloud in terms of its size in the x and y directions. If we assume 
that the cloud has an initial size of zero (treating it as a mathematical point), 
and no initial velocity, we can apply simple kinematics equation that relates the 
speed of the particle with the distance moved: 

v = ~ t (15) 

Now we can replace the velocity, v, and its corresponding differential in equation 
(14) with equation (15) by applying identical method used to obtain equation 
(13), this gives the new distribution: 

Nmx 2 mx 2 
n, = — ea^-^-^cft (16) 
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The equation (16) has the form of a Gaussian distribution, namely: 



1 x 2 

f(x) = —^=exp(-—)dx (17) 
0V27T 2o l 



In the experiment, we can directly measure x which is the variable. We can now 
equate the arguments in the exponentials of the equations (16) and (17): 



2 2 

mx x 



2k B Tt 2 2a 2 



(18) 



Cancelling the x and rearranging for a, we obtain the size of the cloud at a given 
time t and temperature T: 



o = J*^f (19) 
V m 

So in the TOF experiment, we can take a picture of the cloud at different times, 
t, and measure the corresponding cloud size. Then by estimating the gradient of 
the cloud size with respect to time we obtain the gradient which is the argument 
under the square root. This can then be rearranged algebraically to obtain the 
temperature of the cloud. If the cloud has initial size, a , then the time evolution 
of the cloud size is just the sum of the squares: 



°(t) = \J^t* + o* (20) 

Equations (19) and (20) are purely theoretical in form. If we have an optical 
system which images the cloud onto a CCD, then we need to consider the magni- 
fication of the system as well as the pixel size. The cloud size will be measured in 
the number of pixels. Each pixel has a size, l P i xe i, which when multiplied measures 
the actual size of the object produced by the optical system. To get the actual 
size, we then divide by the magnification, M, of the system. Thus, equation (20) 
becomes: 

a actual (t) = lj jf*(t) = lj jf]/^t 2 + a 2 (21) 
2.3 Geometric optics of microscopes 

A typical system that is used in absorption imaging (see Figure 1) consists of 
two lenses that are used to focus an image onto a CCD. Such two lens system 
is generally classified as microscopic as it either magnifies or demagnifies the re- 
sultant image relative to the object. The magnification depends on a number of 
parameters of the system and is not exclusive to the lenses alone. It is therefore 
important to analyse such system and understand its behaviour under various 
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conditions. 



We consider two perfect lenses placed at a distance fi + f% between each other. 
An object is placed at a distance s — fi in front of the first lens (as shown in 
Figure 1). We apply the Thin Lens equation to the first lens: 



1 

Ti 



1 



(22) 



But s — /i, hence we have that s' x = oo. s' : = oo implies that the beam emerging 
on the right hand side of the first lens is collimated. 
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Figure 1: Two lens microscopic system presented in an ideal scenario where the 
object and the image are at the focal lengths of both lenses. 

Now, we want to find the distance S2, from Figure 1 this is just: 

«2 = h + h - s 'i =>■ s 2 = /i + h - oo =^> s 2 = -oo (23) 

As f\ + f'2 00. Now applying the Thin Lens equation to the second lens we 
get: 



11111 

T =-+-^ T =— 

J2 S 2 S 2 J2 



1 



— =>■ h — s 2 

OO s 2 



(24) 



So from this we can conclude that if the object is placed at a focal length of the 
first lens and the distance between the two lenses is the sum of their focal lengths, 
then the final image will form at a focal length of the second lens. 

Conversely, if the beam entering the system is infinity corrected, i.e. s = 00, 
then it will emerge parallel to the optical axis as shown in figure 2. 
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A A + A A 

Figure 2: Ray tracing plot of two lens system. Blue: the object is placed at a 
focal length of the first lens. Red: the beam entering the first lens is parallel to 
the optical axis. 



Note that if the distance between the two lenses is not, /1 + /2, then the emerging 
beam on the right hand side of the OS will no longer be collimated. In other 
words, the beam will be converging or diverging depending on whether the dis- 
tance between, fi + f2, is bigger or smaller. 



The total magnification of an n lens system is just the product of magnifica- 
tions of each lens: 



rriT 



(25) 



Hence, for our two lens system it is just simply: 



rriT 



Si So 
— X — 
Si s 2 



rriT = — 



h 
h 



(26) 



In reality being able to put the object at exactly the focal distance is not nec- 
essarily possible. As a result it is important to understand how the behaviour 
of the system changes when the object is not at the focal point of the first lens. 
Figure 3 shows such situation where the object is placed at a distance, s\ = fi+S, 
where, 5 G V5i. 
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Figure 3: Two lens microscopic system presented in non - ideal scenario where 
the object is not at the focal point of the first lens. 



Once again we apply the Lens Maker's equation on the first lens and rearrange 
to get s x : 



«i = j s 1 = + fi (27) 

si- h o 

Now the object distance s 2 from the second lens is: 

S2 = fl + h ~ s[ ^S 2 = f 2 - f jr (28) 

Applying the Thin Lens equation to the second lens and rearranging for s 2 and 
simplifying it we obtain: 



* 2 = -f(/ 2 -f) ^ 



The final image no longer forms at the focal point of the second lens. However, 
as 5 — > 0, then the equation (28) tends to f 2 as expected. Finally, using equation 
(24) we obtain the total magnification of the system: 

mr ^ x ^M^ x -»(/* ; f> (30) 

si s 2 fi + 5 f2 _L 
The expression above reduces to: 



rriT = —^j- (31) 



From this we can infer that changing the distance between the first lens and the 
object has no effect on the total magnification of the system. 

We now proceed to consider an even more complex case. As in previous part 
we assumed that the object is placed at a distance that is not equal to the focal 
length of the first lens. Now, we will also assume that the distance between the 
two lenses is not f\ + f 2 , but fi + f 2 + a, where, a, is some deviation from the 
ideal distance. Figure 4 depicts the situation. 
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Figure 4: Two lens microscopic system presented in non - ideal scenario where 
the object is not at the focal point of the first lens and the two lenses are not at 
a distance f± + / 2 . 



For the first lens, as previously, the distance s' x is the same, namely: 

' fi 



The new distance s 2 is then: 



S2 = fi + f'2 + a - s[ = fi + f 2 + a 



fi =>- s 2 = f 2 + « 



A 2 



(32) 



(33) 



Applying the Thin Lens equation to the second lens and rearranging to get s 2 we 
get: 



s 2 f 2 (/ 2 + a-?)/2 



S2 - /2 



/1 2 



(34) 



Once again we observe that the new image no longer forms at the focal length of 
the second lends but instead at some distance that depends on the parameters a 
and 5. Once again, in the limits as 5 A a — > 0, the object is formed at the focal 
length of the second lens, as expected. 



The total magnification of the system is: 

/1/2 



ttit = 



a5 - ft 



(35) 



In this case we observe that the magnification is not — /2//1 as in previous case. 
The deviations a and 5 couple. 

The expression above describing the effects of magnification based on variations 
in distance has got some important consequences. First, if the distances a and 
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5 are large, the magnification is unstable for large a5. Hence, if our object's 
position cannot be fixed and can experience some variation in a, then it's re- 
sultant magnification would vary. This instability would have consequences in 
the calculations of temperature of the cloud since our model assumes that the 
magnification is fixed. Therefore, optimising the object to first lens and lens to 
lens distances is of prime importance as it minimises the a5 product and reduces 
the magnification instability even if the position of the object varies from time to 
time in the TOF experiments. In reality, the lens to object distance accuracy is 
harder to achieve, therefore, if the lens - lens distance carries a small deviation, 
5, then the coupling term in equation (35) vanishes as it becomes much smaller 
than ff. 



2.4 Resolution of a two— lens system 

An important quantity in optical design process is the numerical aperture, NA. 
The numerical aperture is a dimensionless quantity that defines the range of 
angles that an optical system can accept or emit light. This is given by: 

NA = nsinO = nsin(arctan(—)) (36) 

Where, n, is the refractive index of the medium, D, is the diameter of the lens 
and, /, is the focal length of the lens. The diffraction limited resolution can be 
estimated using the well known relation: 

x - °^ (37) 

Where, A, is the wavelength of light used. For multiple lens systems, the resolu- 
tion of the system is dictated by the first lens (lens closest to the object that is 
being imaged). The closer the lens is to the object, the higher is the light gath- 
ering power and so the information loss is reduced. In addition, each lens is not 
perfect and therefore has some underlying impulse response or point spread func- 
tion, (PSF), associated with it. The final image for an, n, lens system in Fourier 
plane is given by the convolution between the object and the corresponding point 
spread functions of each lens: 

I(x, y) = 0(x, y) ■ PSF, (x, y) ■ PSF 2 (x, y) ■ . . . ■ PSF n (x, y) (38) 

In equation (35) the order of multiplication matters i.e. it is not commutative if 
the object and the point spread functions are expressed in finite matrix arrays. 



2.5 RIC, OPD and MTF curve characteristics 

There are a number of indicators in computational optics which give information 
on image quality and the overall performance of optical system. Understanding 
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these is therefore key in being able to design an optical system free of defects, 
low aberrations and good resolution. 

Ray intercept curves, or ray trace curves are prime indicators of the aberra- 
tions that are present in the optical design. They can qualitatively tell the type 
of aberrations that exists in the optical design and also give a quantitative rep- 
resentation of how big they are [4]. Figure 5 below shows an example of a RIC 
for a microscope lens system. 
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Figure 5: RIC and OPD curves for an fi = 100mm, fi = 300mm, microscope 
system using OSLO. 

Here the aberrations are 3 rd order spherical aberrations characterised by a cu- 
bic polynomial. The image was placed at an optimum focus that minimises the 
residual blur diameter. Spherical aberrations are characterised by a variation in 
focus position with aperture. 

Optical path difference, (OPD), is another important quantity in optical design. 
It quantitatively describes the difference between a perfect wavefront and a real 
wavefront [4]. Generally, if the OPD is |A, then the system is diffraction lim- 
ited. Here, in the figure above, we can see that the maximum OPD is more than 
3 A, which exceeds the diffraction limit. The distortion of the wavefront varies 
spatially across the optical system. Optical path differences and the distortion 
of the wavefront depend on the quality of the lens and particularly its surface 
smoothness. 

The final and most comprehensive indicator in optical performance that will 
now be discussed is the modulation transfer function, MTF. Figure 6 shows a 
computed MTF curve for fx = 100mm, / 2 = 300m, microscope system using 
OSLO. 
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Figure 6: computed MTF curve for f\ = 100mm, f% = 300mm, microscope sys- 
tem using OSLO. Here the monotonically decaying straight lines represents how 
the MTF would behave in ideal situation and the blue curve represents what the 
performance of the system actually is. 

The modulation transfer function is estimated by comparing the contrast of the 
image with that of the object. Namely, if we have an object and measure its 
minimum and maximum intensity, we can obtain its modulation [4]: 

Modulation = (39) 

I + 1 ■ 

-Lmax ~ x min 

We can then place the object to be imaged by our optical system and perform 
the same procedure of measuring the different intensities and estimate the mod- 
ulation. Finally, to calculate the modulation transfer function we take the ratio 
of the modulations in image and object, namely: 

MTF = M0Dim ° (40) 

By performing such procedure for different size objects, we can compute the en- 
tire spectrum of modulation transfer function based on the object's size. The 
MTF is affected by diffraction, and geometrical aberrations. In addition, the 
MTF curve characteristics can vary greatly across different spatial frequencies, 
i.e. a fairly large object can have a lower modulation than a smaller object of 
higher spatial frequency. 
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The cutoff frequency in the MTF is the highest possible frequency at which 
the modulation goes to zero - the contrast of the system is effectively reduced 
to zero and the different features are of the same intensity, l m in — I-max- This 
allows one to determine the resolution of the system just by reading off the value 
for the cutoff frequency from the MTF, the resolution is given by the inverse of 
cutoff frequency. 

3 Optical design and development 

3.1 Initial stages: establishing constraints and optical per- 
formance requirements 

Before the actual design of the system takes place, there must exist a period 
of development whereby the designer establishes what it wants to achieve with 
the optical system and what sort of constraints limit the optical design. Based 
on these limitations and demands one can then start the process of initial de- 
velopment where mathematics and computer simulations take place to test the 
performance of the system. Finally, given that the system passed all the initial 
development tests it is then sensible to realise the design in a real experimental 
environment and test it. 

The optical system design and optimisation procedure can be summarised in 
the figure 7. We begin to describe each step in the flow chart in more detail [4] . 

1. The first step in the design process should involve acquiring all of the speci- 
fications. These may include optical specifications such as the focal length, 
lens diameter, environmental requirements such as temperature and space 
constraints, performance goal - what resolution or magnification one wants 
to achieve. 

2. The starting point would be a basic configuration that can meet the speci- 
fications of the optical system to some degree. 

3. Establishing variables and constraints. The system variables generally in- 
clude things like radii, thickness, airspaces, glass characteristics, wavelength 
of light used. The lens dimensions are usually restricted by the manufac- 
turer's standards and cannot be varied unless a custom built design is re- 
quested. It is useful to assess the importance and impact of each variable 
in the optical design process and to which extent it can be varied in order 
to achieve the required system specification. 

4. The simple prototype of the system can be programmed using optics soft- 
ware (e.g. OSLO or Zeemax). 
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5. The optimisation process of the system can have a number of requirements. 
For example, the final focus of the image could be optimised such that 
it minimises the optical path differences in the wave front or minimising 
the spherical aberrations by means of introducing lenses with glasses of 
different refractive indices. Having a detailed description of points 1 and 3 
can greatly improve the optimisation procedure. 

6. A number of optical indicators such as the Modulation transfer function 
(MTF), the Ray intercept curves (RIC), the optical path difference curves 
(OPD), rms sport radius and others can provide with a very comprehensive 
overview about the performance of the optical system. It is therefore crucial 
that their properties are well understood. 

7. If the performance of the system is not satisfactory, points 3, 4, and 5 can 
be repeated until the desired performance is achieved. 

8. In reality, it is highly improbable that the system will be built with the 
level of precision and alignment as specified in the simulations. In addi- 
tion, the optical materials and physical conditions may not be perfect and 
homogenous throughout. It is therefore important to investigate the extent 
to which the performance of ones optical system changes with different er- 
rors. If the tolerances appear to be significant, action can be taken early in 
the development process to minimise the risk in propagation of such design 
imperfections in the final stages of assembly. 

9. Provided that the system has passed all the initial tests it can now be 
assembled and investigated further in real life. 




Figure 7: Optical system design and optimisation flowchart. 
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3.2 System constraints in our optical design 

We now proceed to outline some of our optical design constraints and require- 
ments. 

Physical Constraint (Jl. The chamber size put a constraint on the focal length 
of the first lens such that f\ > 137.7mm, this meant that we had to use a mini- 
mum of /i > 150mm (minimum focal length sold by the manufacturer based on 
this constraint - Thorlabs). Of course, a multiple lens system could be designed 
to match more closely to the distance 137.7mm, however, this will certainly add 
several complications in the design process. One of such complications would be 
the care and precision required to assemble multiple lenses systems - this in effect 
would add more degrees of freedom and without the machinery human error is 
an inevitable end. In addition, custom lens designs from the manufacturers is 
very costly and is generally not considered unless the system would have a huge 
improvement on the performance. In short, the best optical system is one which 
satisfies Occam's razor. 

Physical constraint JJ2. The diameter of the imaging beam has also put a 
constraint on the size of our lenses that we had to choose. The imaging beam 
diameter was chosen to be d = 20mm. Hence, to reduce the clipping of the light 
and thus minimise the information loss from the image, the diameter of the lenses 
was chosen to be d[ enses = 50.4mm, this size was limited by the manufacturer. 

Optical constraint (jl. The magnification of the system was chosen to be 
around M > 3. This put a constraint on the choice of our second lens. A lens 
of fi = 500mm was chosen to produce a theoretical magnification of M — 3.33. 
The choice of the focal length of the second lens once again was limited by the 
manufacturers standards (Thorlabs). 

Physical constraint (J3. The above optical constraint has further implicated the 
space constraint. The significant increase in length of the optical system meant 
that it could no longer be built in a linear assembly fashion. Whilst the testing 
could be easily carried out in a linear setup, the assembly and implementation 
of the system in the presence of the UHV chamber could no longer be supported 
by a linear setup. This is fortunately not a significant issue as it can be easily 
overcome by implementing mirror joints which can effectively "fold" the optical 
system into a compact form without suffering degradation in the overall perfor- 
mance. 

There is another alternative which can be used to reduce the effective size of 
our optical system. In the mathematical analysis (Section 2.3) we assumed that 
we are imaging a mathematical point placed at the focal distance from the first 
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lens. This implied that the emerging beam from the second lens was collimated. 
As a result, the second lens could be placed very close to the first one and still 
form an image at a focal point of the second lens. In reality, this can be achieved 
as far as the paraxial approximation is satisfied. However, as discussed in section 
2.3, it adds further complications such as the magnification instability and the 
fact that the laser beam is no longer collimated. 

Optical constraint JJ2. The operating imaging beam wavelength was set to 
A = 780.032nm. This put a constraint on choosing the wavelength operating 
range for our lenses. Thorlabs provides A = 650 — 1050nm operating range. 

Optical constraint (j3. Since our imaging system would be composed of two 
lenses it was important to establish what type of lenses to use. In the initial 
development, piano convex lenses were used. However, it was found that they 
carried large spherical aberrations and optical path differences. The next best 
option was found to be using achromat doublets. These lenses are made out of 
two different types of glass - flint and crown, as a result, they significantly re- 
duced the spherical aberrations and the optical path differences. 

Optical Constraint (J4. Using the information of the known optical constraints 
and the theory on the cloud expansion rates at a given temperature (Refer to 2.2 
Derivation of the temperature of the cloud) we can estimate the constraints which 
our imaging camera must satisfy. Taking the standard mass for Rubidium to be 
m Rb = 85.4678m, where, u = 1.660538921(73) x 10- 27 kg, is the unified atomic 
mass, temperature, T ~ 30fiK, and allowing the time of flight expansion to be 
30ms, we estimate (eqn 19) the size of the cloud, a = 1.61mm. In addition, the 
cloud will fall under gravity for 30ms covering a distance of s = ^gt 2 = 4.41mm. 
Using this information we can estimate the extremum of the cloud edge, this is 
given by the sum of the cloud size and the distance moved: d = 6.02mm. How- 
ever, our system has a magnification that is M = 3.33, and so the resultant image 
projected onto the CCD will be larger by the same factor. Hence, our active sur- 
face area should be no less than 20mm x 20mm. In addition, the pixel size of 
our camera is constrained by the resolution of our imaging system. Therefore, 
to utilise the achieved resolution, one should attempt to choose a camera with a 
pixel size that is smaller than the given resolution, but at the same time not too 
small so that quantum efficiency is not affected too greatly (quantum efficiency 
of photon detection goes down as the pixel size decreases). 

3.3 Development and Simulations 

The initial system development was carried out using OSLO software (Optics 
Software for Layout and Optimization). The lens material specifications and 
dimensions can be readily entered into OSLO and optimised using inbuilt algo- 



16 



rithms. 



The simulation procedure involved the following steps: 

1. Build the first lens. Optimise the lens distance such that the distance be- 
tween the object and the lens produces a collimated beam on the right hand 
side. The optimisation procedure would be carried out using aberration 
operands, more specifically, using, PU, or paraxial axial ray optimisation 
operand. We have included the chamber window in the optical simulation 
to account for any effects due to the extra optical components. The beam 
radius was set to, r = 10mm, wavelength, A = 780.032nm. The achromat 
doublet was placed with the flatter end (large radius of curvature) to face 
the object. This reduces the level of spherical aberrations and distortion of 
the wavefront. 
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Figure 8. Achromat doublet with the mathematical point at the focal point. 
The optimisation was carried out using paraxial axial ray optimisation. 

2. A second lens is placed at an arbitrary distance from the first lens. The 
object is now placed at an infinity (the rays are collimated). Once again, 
we use Paraxial-axial ray optimisation algorithm to optimise the distance 
fi + fi so that the emerging beam is collimated on the right hand side of 
the optical system. The second lens is placed such that it faces the first lens 
with a surface of smaller radius of curvature. Figure 9 displays the setup 
in OSLO. 
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Figure 9: A pair of achromat doublets with image placed at infinity. Here 
fi = 150mm, and /2 = 500mm. 

From this simulation we can observe that the emerging beam from the 
first lens has a diameter that is larger than the second lens. This implies 
that there will be some clipping of light involved. Clipping of light results 
in the loss of information about the image and a lower signal to noise ratio. 
To remedy this there are a number of options. First, we could increase the 
diameter of the second lens so that all of the light is collected. However, in 
our case we cannot do that since the manufacturer does not offer lenses with 
a bigger diameter. Secondly, we could sacrifice the magnification stability 
by reducing the distance between the two lenses. However, in our case that 
could seriously affect the calculations of the atom number and temperature 
of the cloud. 

3. Finally, the optimised distances are combined in a system were the object is 
placed at a focus of the first lens and the distances between the two lenses 
are placed at the optimised distance fi+fe. The final image is focused using 
an Autofocus optimisation option for minimum optical path difference for 
monochromatic light. Figure 10 shows the final optimised setup and the 
table below summarises the optical component distances. 
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Figure 10. Optimised achromat doublet microscopic system for ultracold 
atom imaging. 



Table 1: Simulation results 



s t = fx, (mm) 


fx + h, ( mm ) 


4 = fo, (mm) 


146.5 


654.8 


494.3 
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We now begin to analyse the performance of our optical system by looking at the 
parameters such as the RIC, OPD, MTF. Figure 11 displays the RIC and OPD 
curves. 

From the OPD curve we can see that the wavefront distortion is nearly — |A 
- which is close to the Rayleigh limit. In addition, the spherical aberrations are 
3 rd order spherical aberrations. If we were to compare this level of aberrations 
with an identical achromat doublet system, but with their faces reversed (i.e. 
smaller radius of curvature facing the object), we would find that the aberrations 
are 6 times smaller with an OPD maximum being 4 times smaller in the first 
system compared to the second one. 
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Figure 11: Optical Path difference and Ray intercept curves for the achromat 
doublet system, fx = 150mm, f'2 = 500mm. 

The modulation transfer function of our system has the following behaviour de- 
picted in figure 12. The MTF curve of our system follows closely the behaviour of 
the ideal system which indicates small geometric aberrations. This is confirmed 
by the OPD and RIC curves. 

The resolution of the system can be estimated by reading off the cutoff frequency 
from the MTF curve, namely: 

1 

^cutoff 

Taking the cutoff frequency to be approximately equal to v C utoff — 70 cycles 
mm~ l . With this value, we can estimate the theoretical resolution over one cycle 
which is defined as a distance between the brightest point and the darkest point. 
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Figure 12: MTF curve for the achromat doublet system, fx = 150mm, f 2 = 
500mm. Here the black curve depicts the MTF of an ideal lens system and the 
blue curve represents how our lens system should behave in reality. 

This gives us a theoretical estimate for the resolution of the system: 

11 

x m in = = ^ = 14.3/im (42) 

"cutoff 'U 

If we were to base our resolution estimate on equation (37), then we would find 
that it yields: 

Xf 780.032nm • 150mm . , 

x min = 1.22-^- = 1.22 = 4.76um 43 

D 30mm p v ; 

In our case, the best resolution that can be achieved is limited by the chamber 
window (chamber window clear aperture diameter D = 30mm). 

We can observe that the two estimates show a substantial difference. This can be 
explained, at least in some part, by the fact that the spatial resolution estimate 
using equation (43) does not take into account any geometrical aberrations that 
may be present in the system. 

The magnification of the system given by OSLO was M = 3.31. 



3.4 Tolerances and Errors 

Whilst it is possible to achieve a perfectly optimised system in the simulation, in 
real experimental conditions this is not the case. The optimum distances between 
the optical components and object/image may not necessarily be achieved. As 



20 



a result, it is important to understand how well our optical design can tolerate 
such imperfections or errors and how much it can affect the image quality. 



First we proceed to investigate the tolerancing effects of the object being not 
at the optimum focus of the first lens. Increasing the object-lens distance we 
observed astigmatism in the ray intercept curves. The wavefront distortion in- 
creased significantly with the wavefronts having largest optical path differences 
near the edges. The performance described by the MTF curve also showed nega- 
tive effects, though the resolution of the system remained unchanged. Figure 14 
summarises the tolerancing effects. 
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Figure 13: System performance when we place the object at s = fi + S 
150 + 5mm. 



Next, we move the object closer to the first lens. This time, we observe astig- 
matism and defocus in the ray trace curve. The wavefront distortion follows the 
same profile but is reflected in the x-axis. The MTF curve once again shows a 
decrease in performance, though no noticable change in resolution was observed. 
Figure 14 shows the results. 
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In both cases the system displays similar characteristics, however, the MTF curves 
differ due to different profiles of geometric aberrations. We now proceed to vary 
the distance between the two lenses as well as the distance between the object 
and the lens. 
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Figure 14: System performance when we place the object at s = /i + 8 = 
150 — 5mm. 



The tolerancing of the system when the lens-object and lens-lens distances were 
varied simultaneously did not show any noticeable changes in geometric aberra- 
tions, i.e. the geometric aberrations and the performance of the system remained 
identical to the case were only object-lens distance was varied (in each case we 
varied the distances by up to ±5mm). However, a change in magnification was 
observed. The magnification varied by few parts in hundreds of decimal point. 
This is expected from the analysis carried out in section 2.3. 

Now we look at how the system performance varies when the lens-camera dis- 
tances are not at their optimum points. Once again these are varied by ±5mm. 
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Figure 15 depicts the results. 



We can observe that the aberrations and wavefront distortion is much smaller 
relative to the errors produced by the first lens. The aberrations in the sys- 
tem remain 3 rd order spherical aberrations. The optical path differences have 
increased, though are by a factor of 10 smaller than the ones produced by the 
first lens. The MTF curve seems to have poor performance in the intermediate 
regime, however, for small features it performs as an ideal system. 
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Figure 15: System performance when we place the camera at f2 + $ = 500 + 5mm 

Qualitatively, when the distance between the lens and the camera is decreased 
by the same amount, the aberrations remain spherical and the OPD maximum 
is reduced by a factor of 2. The MTF curve performed poorly across different 
spatial frequencies. 



When both object-lens and lens-camera distances were varied the aberrations 
and optical path differences including the MTF performance are completely dom- 
inated by the effects produced by the first lens. 
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4 Building and Testing 



4.1 Lateral Imaging 

The first part of imaging system testing was in the lateral imaging configuration. 
This is where the collimated laser beam incident on the object is completely 
parallel to the optical axis. We have used single mode optical fibre fed into a 
collimator as an imaging beam. 

The setup configuration was realised on an optical bench with optical components 
being mounted on rails and linear translation stages. As an imaging target, we 
have used a 1951 USAF target. For imaging, we have selected a Guppy Pro F- 
125 camera with 1292 x 964 resolution and 3.75/im pixel size. Figure 16 depicts 
the experimental setup. 



/i = ISO ± 5mm f 2 = 500 + 5mm 




/ * > 

/ Fin» sdjuslmsnl lif*ar _ Fine jdjustmsnf linear 

/ translation stagr (grey) fx + ft = 650 + 5WW Tr*nl1*(l0n *<*gd [K'Cy) 



Rail: 

Figure 16: Experimental setup for testing the lateral imaging system (not to 
scale). We used 7T175-50 - Aluminium Translation Stages with 50 mm travel 
range. 

With this experimental arrangement we have obtained the following image shown 
in figure 17. The image has some interference effects caused by coherent and 
monochromatic laser light. Some of the interference is produced by small debri 
on lenses and the target. 

To determine the magnification of the system, we selected a square feature 
(marked blue on figure 17) which was the 2 nd element in group 4. The square 
covered 2.5 cycles, and it had a spatial frequency of 17.9 cycles mm' 1 (as read 
from the standard 1951 USAF size chart). From this it was possible to calculate 
the size of the square over 2.5 cycles which gave 0.140mm. The 0.140mm is the 
size of the square if the magnification of the system is M = 1. Using Image J 
software the length of the square was measured to be 126 ± 1 pixels in x and y 
direction (to nearest pixel). 
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Figure 17: 1951 USAF target using the lateral imaging configuration. The red 
square indicates the smallest feature that could be resolved on the target whilst 
the blue feature was used to estimate the magnification of the system. 



Knowing the pixel size, (3.75/im), the size of the magnified image could be cal- 
culated. It is simply the product of pixel size and number of pixels measured 
namely: length = 126 ± Ipixels x 3.75fim = 0.473 ± 0.004mm. The error has 
been obtained by taking repeated measurements of length and estimating the 
standard error on the mean. This error calculation is not extensive as it does not 
take into account any systematic, rounding or alignment errors. Now, to estimate 
the magnification of the system, we divide the experimentally estimated size by 
the theoretical estimate of magnification 1: 



M 



length 



exp 



length 



theo 



0.473mm ± 0.004 
0.140mm 



3.38 ±0.03 



(44) 



Now the smallest feature that could be resolved in the x and y direction (marked 
red) was element 5 in group 6. From the standard 1951 USAF target this had 
a spatial frequency of 102 cycles mm -1 . With this value, the resolution of the 
system is: 



1 



1 



^cutoff 102cyclesmm 1 



9.80/im 



(45) 
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The standard convention in estimating the resolution using a 1951 USAF target 
requires the features in horizontal and vertical direction to be resolved and the 
resolution is defined over one full cycle which is the distance between the bright 
and dark fringe. 

The experimental estimates of resolution and magnification are relatively sim- 
ilair to the estimates obtained in the theoretical calculations and simulations, 
though there are some differences. 

4.2 Grazing Incidence Imaging 

In grazing incidence imaging, the collimated imaging beam is reflected off an 
atom chip (in our case mirror). If an object is close enough to the mirror, a vir- 
tual image is produced and the final image consists of two objects. The distance 
between the two objects allows one to determine the chip-object distance. 

To simulate an atom chip we have used a 50.8 mm protected gold mirror from 
Thorlabs. It was mounted parallel on a separate rail. For fine adjustment, the 
mirror was attached to a 6mm range, 7T67-6 - Stainless Steel Translation Stage 
from Standa. 

It was not practical to implement the 1951 USAF target since the smaller features 
were far away from the mirror. Instead, the target was replaced by a single hair 
of 0.05 ± 0.01mm thickness. The imaging beam was tilted using fine adjustment 
screws. A protractor was used to measure the angle of inclination. Figure 18 
depicts the experimental setup. 
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Figure 18: Experimental setup to test grazing incidence imaging (not to scale). 

In grazing incidence imaging, we tilted the beam by 1° — 2° degrees. The re- 
sultant images obtained are depicted in figure 19. 
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Figure 19: Grazing incidence imaging. Left image, 6 = 1° ± 0.5°. Right image, 
6 = 2° ± 0.5°. The interference fringes in the images are produced due to diffrac- 
tion of light around the edges of the mirror. 

It must be noted that the mirror blocks some of the light incident on the ob- 
ject and so the signal to noise ratio increases. 

To determine which one of the images is virtual and which one is real, the object 
was moved and the corresponding motion of the two images observed. The real 
image moves in the same direction as the object moved, whereas the virtual im- 
age moves in the opposite direction. 

The distance, d, between the virtual and real image can be measured and used to 
calculate the distance of the object away from the chip. In terms of geometrical 
optics, some of the imaging beam is first incident on the chip and then reflected 
back onto the image which produces a shadow - this is the real image. In the 
second case, the imaging beam is first incident on the object which then casts a 
shadow on the chip which is then reflected off it generating a virtual image. 




Figure 20: Ray propagation in grazing incidence imaging [1]. 
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Using trigonometry, the distance, d, can be estimated using [1]: 

d = 2hcos9 (46) 

Hence, to find the chip-object distance, h, we can just rearrange the equation 
above. It must be noted, that if the object gets too close to the chip, then the 
virtual and real images will merge together. This of course depends on the size 
of the object being imaged as well as the resolution of the system. 

5 Conclusion 

In this report I have discussed the technical details on how to develop and build 
an absorption imaging system for ultracold atoms. The optical design, whilst rel- 
atively simple, has numerous capabilities in probing and extracting information 
from cold atoms ranging from the temperature calculations to the distances from 
the chip. However, this simplicity appears to be superficial and various technical 
intricacies begin to emerge when system constraints are taken into account. 

It is evident that, in order to improve the performance of the optical design 
one should attempt to ease the optical and physical constraints as these pri- 
marily limit the choice of the optical components and ultimately determine the 
resolution of the system. 
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